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ABSTRACT 

Unli ke Trojans, horseshoe coorbitals are not generally con sidered to be long-term sta- 
ble permott fc Murravl Il981al iMurrav fe Dermottl Il999l) . As the lifetime of Earth's 
and Venus's horseshoe coorbitals is expected to be about a Gyr, we investigated the 
possible contribution of late-escaping inner planet coorbitals to the lunar Late Heavy 
Bombardment. Contrary to analytical estimates, we do not find many horseshoe ob- 
jects escaping after the first f 00 Myr. In order to understand this behaviour, we ran 
a second set of simulations featuring idealized planets on circular orbits with a range 
of masses. We find that horseshoe coorbitals are generally long lived (and poten- 
tially stable) for systems with primary-to-s econdary mass ratios larger than about 
1200. This is consistent with the results of Laughlin & Chambers) (|2002f ) for equal- 
mass pairs of coorbital p lanets and the instability of Jupiter's horseshoe companions 
(jStacev fc Connors] 120081) . Horseshoe orbits at smaller mass ratios are unstable be- 
cause they must approach within 5 Hill radii of the secondary. In contrast, tadpole 
orbits are more robust and can remain stable even when approaching within 4 Hill 
radii of the secondary. 

Key words: celestial mechanics - planets and satellites : dynamical evolution and 
stability - methods : numerical - Earth -planets and satellites : individual : Venus. 



INTRODUCTION 



A horseshoe orbit is a type of coorbital motion where two 
bodies maintain the same average distance from the central 
object, while librating around points 180° apart in longi- 
tude. If one body is much more massive than the other, 
the smaller object's trajectory resembles a horseshoe when 
plotted in the frame rotating with the bodies' average mean 
motion. Horses hoe orbits were a li ttle-known theoretical pre- 
diction (|Brownlll91ll ; lRabdll96lr ) until it was realized that 
Saturn 's moons Janus an d Epimetheus actually follow such 
paths (|Smith et alJll980h . 

The most detailed theo retical exploration o f horsesh oe 
orbits was undertaken by iDermott fc Murravl (|l981al lbT). 



approach, IDermott fc Murravl l|l981al ) derived the following 
formula for the horseshoe lifetime 



5/3 



(1) 



where T the orbital period and \i is the secondary-to- 
primary mass ratio. While this formula was originally de- 
rived as a first-order estimate of how long the horseshoe mo- 
tion may last (S. Dermott, 2010, personal communication), 
it was later presented as a relatively firm upper limit on 
the l if etime of such configuratio n s (e.g. IMurrav fc Dermotd 
1 19991 ; IStacev fc Connors! 12008) ; IChristou fc Ashed 120111 ) 
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Using an analytical approach, 
could not confirm that succes- 
sive horseshoe loops repeat exactly, even in the restricted 
problem. Instead, they derived an estimate of possible drift 
during one horseshoe cycle, and extrapolated the lifetime of 
the horseshoe configuration assuming that these drifts over 
long periods of time behave like a random walk. With this 
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Note that this relation indicates that the observed horse- 
shoe coorbitals Janus and Epimetheus can last longer than 
10 Gyr, while Jovian horseshoes should be unstable on Myr 
times cales (as confirmed numerically by IStacev fc Connors] 
2008). Given that this relation is based on a relatively simple 
extrapolation, the ultimate intrinsic stability of horseshoe 
orbits must be explored using direct integrations which were 
not possible in the Voyager era. 

The initial motivation for this study was the question 
of the source of late luna r impactors, usually termed the 
Late Heavy Bombardment (|Chapman. Cohen, fc Grinspoonl 
2007). The formation of the Moon's giant Imbrium basin 
3.85 Gyr ago, accompanied by many other smaller im- 
pacts, indicates that a population of impactors was present 
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Figure 1. Remaining stable Trojans (squares) and horseshoe coorbitals (circles) of Earth after (a) 10 Myr, (b) 100 Myr, (c) 300 Myr 
and (d) 700 Myr. Pluses indicate unstable particles. At t=0, all orbits had e=0.05 and the same longitude of the node and pericenter as 
Earth. 



in near-Earth space about 700 Myr after the formation 
of the planets. While non-coorbital inner Solar System 
planet esimals should ha ve been mostly cleared out by 
then dBottke et al.l 120071 ') , Eq. [T] indicates that horseshoe 
coorbitals may linger for a Gyr or so. Eq. [1] gives the 
horseshoe lifetimes of 1.6 Gyr for Earth and 0.7 Gyr 
for Venus, the right order of magnitude to explain the 
Late Heavy Bom bardment (and consistent with previous 
numerical results; iTabachnik fc Evans! [2000). Furthermore. 
IScholl. Marzari, fc Tricaricd (120051 ) found that Venus Tro- 
jan£| may be also unstable on Gyr timescales. While the 
question of late impactors on the Moon has motivated our 
original numerical experiments (Section 2), our study was 
eventually broadened to address the intrinsic stability of 
horseshoe coorbitals in general (Section 3). 



2 COORBITALS OF EARTH AND VENUS 

As we were initially interested in both horseshoe and tadpole 
coorbitals, we decided to place our test particles 60° ahead 
of the planet in longitude, consistent with both types of tra- 
jectories. One thousand particles were spaced on a uniform 



1 Here we use "Trojan" to mean a tadpole coorbital of any celes- 
tial body. 



grid, with 25 different starting inclinations and 40 differ- 
ent initial semimajor axes. For both Earth and Venus coor- 
bitals, the semimajor axis range was 0.01 AU, starting at the 
planet's mean distance, and the inclination range was — 36° 
(measured from the planet's orbital plane, with the same 
node). This range of initial distances covers about one Hill 
radius from the planet's own orbit, which is the expected size 
of the coorbital region (see Section 3) . All particles had ini- 
tial eccentricities of 0.05 and the same longitude of pericen- 
tre as their parent planet. The eight planets were started on 
their initial positions for January 1st, 2000, and their initial 
vectors and masses were obtained from JPL's HORIZONS 
online service. To integrate the particles' orbits, we used the 
Swift-RMVS4 integrator (provided by Hal Levison), based 
on pr evious successful versions of Swift (|Levison fc Duncan! 
Il994h. Swift- RMVS4 is a "dem ocratic" symplectic integrator 
(cf. IWisdom fc Holman1ll99ir) that uses canonical heliocen- 
tric coordinates and is capable of resolving close encounters 
between test particles and planets, and the RMVS4 vari- 
ant is notable for planets being insensitive to the presence 
of test particles (in previous versions, systems having iden- 
tical planets but different massless particles would slowly 
diverge). With RMVS4, the whole grid would experience 
the same planetary perturbations, even though the compu- 
tation is divided between 40 different processors. We used 
a time step of 3 days for all runs. The Earth and Venus 
grids were integrated for 700 Myr and 1 Gyr, respectively. 
Particles that collided with the Sun or a planet, or reached 
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Figure 2. Remaining stable Trojans (squares) and horseshoe coorbitals (circles) of Venus after (a) 10 Myr, (b) 100 Myr, (c) 300 Myr 
and (d) 700 Myr. Pluses indicate unstable particles. At t=0, all orbits had e=0.05 and the same longitude of the node and pericenter as 
Venus. This integration was continued to lGyr with little change from panel (d). 



100 AU from the Sun, were removed from the simulation. 
All computations in this section were completed on Harvard 
University's computing cluster " Odyssey" . 

Fig. [T] shows the results of our integrations of Earth 
coorbitals. Panels a-d show the distribution of surviving 
horseshoe and tadpole objects on the grid of initial condi- 
tions, at 10, 100, 300 and 700 Myr, respectively. Fig.[2]shows 
the results for the Venus coorbital grid, in the same format. 
The inner part of both planets' coorbital region is taken 
up by Trojan librators, while the oute r part is populated 
by h orseshoe coorbitals, as expected (|Murrav fc DermottJ 
1999)0 In both cases, Trojans appear to be unstable 
at specific inclination s , an effect of secular resonance s 
ijBrasser fc Lehtol 1200 j: IScholl. Marzari. fc Tricaricd l2005n . 



As found before by IScholl. Marzari. fc Tricaricol l|2005l ) for 
Venus, horseshoe coorbitals are most stable at low inclina- 
tions. Contrary to our expectations from Eq. [I] many horse- 
shoe coorbitals appear to be stable at longer timescales, and 
there is little change in the number of stable horseshoe ob- 
jects between 300 Myr (panels c) and 700 Myr (panel d) 
in both figures. The evolution of the two grids is summa- 
rized in Fig. [3] which shows the number of survivors in each 



2 Particles found over more than 180° of mean longitude in the 
frame rotating with the planet's mean motion were considered to 
be on horseshoe orbits, while those confined to a smaller range of 
longitudes were assumed to be Trojans. 



population over time. The number of surviving horseshoe 
coorbitals of both planets seems to plateau after a couple 
hundred Myr, and there are barely any escaping horseshoe 
particles after 500 Myr. Trojans appear to have more of a 
long-lived but unstable tail, but also appear incapable of 
adding sufficiently to the Late Heavy Bombardment. 

It is interesting to compare our results to the orbit of 
the most long- lived temporary horses hoe coorbital of Earth, 
2010 S016 (jChristou fc Asherl [ioil] ). It has a somewhat 
higher eccentricity that our grid (0.075 vs. 0.05), and its 
semimajor axis semi-amplitude (0.004 AU) and inclination 
(14.5°) place it on the edge of the peninsula of instability 
near the centre of Fig. [T^, which becomes a vertical gap in 
panels b-d. Therefore, our integrations are consistent with a 
relative ly long-lived, but ultima tely unstable nature of this 
object |Christou fc Asherl lioTll ). 

Our results diver ge from those of 

IScholl. Marzari. fc Tricaricd (|2005l ) for Venus Trojans. 
They identified a subset of thirty most stable low- 
inclination Trojans of Venus and followed them for 1 Gyr, 
during which time all of the particles became unstable. 
Figs. [2] and [3] show that there is still a significant number 
of surviving Ve nus Trojans at 1 Gyr, which were unlikely 
to be missed bv lScholl. Marzari. fc Tricaricd (|2005l ). Apart 
from the timescale, we also find different profile of Trojan 
loss for Venus: our results show a slow exponential-like 
dec ay, while the population dr o ps ap proximately linearly 
for IScholl. Marzari. fc Tricaricd (|2005l ). Despite the small 
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Figure 3. The number of surviving particles on tadpole (squares) 
and horseshoe (circles) orbits as a function of time, for the coor- 
bitals of Earth (top) and Venus (bottom). 



numbers of particles, we find that our results are signif- 



icantly different statistically. IScholl. Marzari. fe Tricaricol 
(2005) also used Swift, so it is unlikely that the difference 
is an artifact of our integration schemes. 

We hypothesize that there is a hidden physical parame- 
ter that makes our Venus Tr o jans m ore stable than those of 
IScholl. Marzari, fc Tricaricol |2005). We explored the possi- 
bility that th e diver gence of planetary orbits due to chaos 
i|Laskail lliffli , | 2008}) may lead to different incarnations of 
the Solar System in our simulations, which could have conse- 
quences for stability of coorbitals. In particular, the strength 
of the normal mode associated with the eccentricity of Mer- 
cury can vary significantly over time, possibly also affecting 
the secular behavior of Venus Trojans. In order to test this, 
we integrated a subset of relatively stable coorbitals of both 
Earth and Venus ten times with slightly different initial con- 
ditions (achieved by changing the distance of Earth from 
the Sun by a very small amount). While the integrations 
did diverge, there was no real difference in the stability of 
coorbitals over 500 Myr. We therefore think that the secular 
state of the system is not the reason for the discrepanc y be- 
tween the results of lScholl. Marzari. fc Tricaricd |2005l ) and 
our own. One remaining possibility is that planetary mean- 
m otions may be slightly different betw een our runs and those 
of lScholl, Marzari. fc Tricaricd |2005l ) . We did not take care 
to avoid "cold start" issues and match the mean motions 
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Figure 4. Initial radial displacement Aa (in planet's Hill radii) 
of stable (circle) and unstable (pluses) particles as a function of 
planet mass in our idealized one-planet, 10 8 -orbit simulations. 
For each planet mass, twenty particles were evenly distributed 
between the Trojan- horseshoe boundary and one Hill radius be- 
yond the planet's orbital distance, 60° ahead of the planet in 
longitude. Most (but not all) of these orbits are horseshoes. 



of planets precisely to the observed values (as symplectic 
integrators take large steps, correct initial conditions can 
generate an orbit that is slightly different from the correct 
one). Therefore, if some combination of planetary motions is 
near- resonant with Trojan librations (cf. iNesvornv fc Donesl 
2002), it would be very sensitive to small changes in the plan- 
ets' orbital periods, and may have very different effects in 
our simulations and those of (|Scholl. Marzari. fc Tricaricd 
2005). However, this issue is too tangential to explore here 
in more detail, as we are mainly concerned with the stability 
of horseshoe coorbitals. We therefore leave this discrepancy 
unresolved. 



3 COORBITALS IN SINGLE-PLANET 
SYSTEMS 

While the coorbitals may not be a source of late lunar im- 
pactors, we are still interested in the issue of their intrinsic 
stability. Some horseshoe coorbitals of Earth and Venus ap- 
pear to be long-term stable, in defiance of Eq.[T] In order to 
study their intrinsic stability, we decided to run a number 
of numerical experiments in planar systems with only one 
planet on a circular orbit. By eliminating other planets we 
can make sure that any instabilities are not due to their per- 
turbations, and also significantly speed up our integrations. 
First we ran seven simulations in which our lone planet had 
mass of v = io- 2+fex0 - 5 , where k = 1 - 6. Twenty test par- 
ticles were placed 60° ahead of the planet, uniformly spaced 
between the Trojan-horseshoe boundary at Aoi = «o\/8/3/i 
and the Hill radius A02 = ao(/^/3) 1//3 (the coorbital region 
is actually about 20% wider than a Hill radius, but here we 
are only interested in the most stable particles). The planet 
was assumed to orbit a solar-mass star at 1 AU, and the 
timestep was 10 days. All particles were followed for 10 s 
orbits. Results are plotted in Fig. U 

The boundary between stable and unstable particles 
is not clear-cut in either initial semimajor axis or planet 
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Figure 5. Smallest distance in mean longitude from the planet 
for stable particles with initial conditions plotted in Fig 4. The 
particles were integrated for 10 s orbits. The theoretical boundary 
between tadpole and horseshoe orbits at closest approach (23.5°) 
is plotted by the dashed line. Most (but not all) particles above 
the line are on tadpole orbits. 



mass. As the mass range fi = 10~ 4 — 10~ 3 was found 
to contain the stability-instability transition, we filled in 
that region with more closely spaced integrations. Note that 
some of the innermost stable particles along the horseshoe- 
Trojan boundary actually follow tadpole orbits. In any 
case, it appears that all horseshoe coorbitals become un- 
stable once the planet mass reached about one Jupiter mass 
(/i ~ 0.001). After we o btain ed this result, we realized that 
lLaughlin fc Chambers! |2002) reached a similar conclusion 
when they studied horseshoe configurations of two equal- 
mass planets. They found that two Saturn-sized planets can 
make a stable pair (at /i = 3x 10 -4 each), but more massive 
ones cannot. It appears that the stability boundary depends 
only on the total mass in two bodies, regardless of how it 
is distributed. This is, in fact, expected from the equations 
of motion of the system which pr edict dependence on only 
the sum of the coorbital masses (|Yoder et al.lll983l ). How- 
ever, given how "fuzzy" the boundary appears in Fig. [4] it 
is very likely that some other parameter of the horseshoe 
contributes to its stability. 

In Fig. [S] we plot the minimum approach angle of sta- 
ble particles' guiding centres against the planet's mass[f| It is 
now clear that stable horseshoe particles' guiding centres ap- 
pear to be restricted to a relatively well-defined area which 
disappears for fi > 10~ 3 . The upper boundary of this re- 
gion is almost independent of planet m ass, and corresponds 
to the Trojan-horseshoe dividing line. iMurrav fc DermottJ 
(1999) show how, when crossing the planet's heliocentric 
distance at the closest approach to the planet, the guiding 
centres of horseshoe particles must always pass closer than 
23.5° of the planet, while the Trojans stay further away. Our 
integrations agree with this result, as most stable particles 
above the dashed line in Fig.[5]are actually on tadpole orbits 



The guiding centre is a point that has the same semimajor axis 
and mean longitude as the particle itself, but with zero orbital 
eccentricity. In a rotating reference frame, the guiding centre is 
the centre of the particle's epicycle. 



(they were initially placed exactly on the Trojan-horseshoe 
boundary). 

The bottom boundary of the stable region is curved, 
and appears to be a relatively weak function of planet mass. 
As planetary encounters usually scale with the Hill radius 
of the planet, we re-plotted the same data in Fig. [6] where 
the minimum distance between the particle and the planet 
was expressed in Hill radii (Rh = ao^/S) 1 / 3 ), rather than 
degrees. The minimum distance for stable horseshoe coor- 
bitals appears to have a lower boundary at 5 Rh, which 
keeps their guiding centres at least 5.5 Rh from the planet. 
This seems to be the determining factor for the lack of 
stability for horseshoes of Jupiter-mass and larger planets: 
even the innermost horseshoe trajectory for planets of this 
size experiences a relatively close approach to the planet at 
the horns of the horseshoe. This critical mass for horseshoe 
stability (/i ~ 1/1200) is the real expl anation for the in- 
stabi lity of Jupiter horseshoe coorbitals l|Stacev fc Connors] 
2008), and the stability of low-mass systems like Janus and 
Epimetheus, rather than E q. [1] Furthermore, as noted by 
(IDermott fc Murravlll981ah and demonstrated in Fig. 3] a 
larger part of the coorbital region is taken by horseshoes for 
smaller secondary mass. 

This 5 Rh (or 5.5 Rh for guiding centres) minimum 
close-approach distance applies only to horseshoe coorbitals. 
We also integrated a number of Trojan orbits for planets of 
larger mass (squares Fig. [5} and found that Trojans of high- 
mass planets can approach within 4 Rh of the planet and 
still be stable over 10 8 orbits. It is not clear if the Trojans are 
intrinsically more stable than horseshoes when subjected to 
the same planetary perturbations, or the minimum distance 
simply decreases with increasing planet mass. Not coinci- 
dentally, the stable region for Trojans in Fig. |S] disappears 
at about the same mass rati o at which the Lagrange point 
becomes unstable, \x = 0.04 (jMurrav fc Dermottll 19991 ). 

In order to check that our results are not integrator- 
dependent, we ran a second set of simulations with some 
of the same initial con dition propagated using the HN Body 
symplectic integrator |Rauch fc Hamiltor]|2002l . [20T2T ). Fig. 
[7] shows the results of these integrations, which also ran for 
10 8 orbits, but the closest approach distance was calculated 
every timestep (in Fig. [6] the minimum distance was only 
computed for output times 0.1 Myr apart). While some sta- 
ble horseshoe coorbitals do occasionally approach as close 
as 4 Rh to the planet, the region of contiguous horseshoe 
stability ends at 5 Rh, in agreement with Fig. [5] 



4 CONCLUSIONS 

Using direct numerical integrations, we find that a substan- 
tial number of horseshoe coorbitals of Earth and Venus ap- 
pear to be long-term stable. Some of the low-inclination Tro- 
jans of these two planets are also stable over 1 Gyr, but the 
escapes from the stable region continue throughout the sim- 
ulations, likely due to secular interactions with the planet s 
IIBrasser fc Lehtol |2002| ; IScholl. Marzari. fc Tricaricol l2005h . 
More work is needed to ascertain if any of the Trojans of 
Earth and Venus are stable over the Solar System's lifetime, 
and if not, what are the timescales and exact mechanisms of 
their instability. In any case, the lack of observed primordial 
coorbitals of Earth or Venus, together with our finding of 
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Figure 6. Smallest distance from the planet (in planetary Hill 
radii) for the stable particles' guiding centres (top panel) and the 
particles themselves (bottom panel). Circles plot stable particles 
from Figs. 4 and 5 and squares plot a separate set of stable Trojans 
of higher-mass planets. Dashed and solid lines plot the distance 
of the tadpole-horseshoe boundary and the triangular Lagrange 
point, 23.5° and 60° from the planet in longitude, respectively. 
No stable horseshoe coorbitals approach within 5 Hill radii of the 
planet, while Trojans of massive planets can approach within 4 
Hill radii. 



their relative stability, indicate that these objects were not 
an important contribution to late bombardment on Earth 
and the Moon. 

We also integrated a number of horseshoe coorbitals 
in planar, circular one-planet systems. We find that only 
planets smaller than about one Jupiter mass (p = 1/1200 
or smaller) can have stable horseshoe coorbitals. The main 
reason for this stability boundary is that the horns of the 
horseshoe approach too close to the planet (in terms of the 
planet's Hill radii) for more massive planets. This result i s 
in agreement with the work of lLaughlin fc Chambers! (|2002T ) 
on equal-mass planet pairs, the known instability of Jupiter 
horseshoes, and the long-term stability of the only observed 
horseshoe pair, Janus and Epimetheus. 
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